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Qh! Abstract 

^^ , We propose a new method for constructing exact solutions to nonlinear delay reaction-diffusion 

0^ ' equations of the form 

ut = kuxx + F{u,w), 

where u — u{x, t), w = u{x, t — r), and t is the delay time. The method is based on searching for so- 
^^ I lutions in the form u ~ J2n=i Cn {^)Vn (t), where the functions ^„ (x) and r/n {t) are determined from 



a 



additional functional constraints (which are difference or functional equations) and the original delay 
partial differential equation. All of the equations considered contain one or two arbitrary functions of 
a single argument. We describe a considerable number of new exact generalized separable solutions 
and a few more complex solutions representing a nonlinear superposition of generalized separable 
and traveling wave solutions. All solutions involve free parameters (in some cases, infinitely many 
<^ I parameters) and so can be suitable for solving certain problems and testing approximate analytical 

(^ . and numerical methods for nonlinear delay PDEs. The results are extended to a wide class of non- 

[■■"-■ ' linear partial differential-difference equations involving arbitrary linear differential operators of any 

■n[ , order with respect to the independent variables x and t (in particular, this class includes the nonlin- 

ear delay Klein-Gordon equation) as well as to some partial functional differential equations with 



vn 



■^ . time-varying delay. 
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vh '• 1 Introduction 

Nonlinear delay partial differential equations and systems of coupled equations arise in biology, bio- 
physics, biochemistry, chemistry, medicine, control, climate model theory, ecology, economics, and 
many other areas (e.g., see the studies llll-lllll and references in them). It is noteworthy that similar equa- 



tions occur in the mathematical theory of artificial neural networks, whose results are used for signal and 



image processing as well as in image recognition problems 1112142 111 . 



The present paper deals with nonlinear delay reaction-diffusion equations llL |3l [nl |22|] of the form 

Ut = kuxx + F{u,w), w = u{x,t — T). (1) 

A number of exact solutions to the heat equation with a nonlinear source, which is a special case 



of equation ([T]) without delay and with F{u,w) = f{u), are listed, for example, in ll23M29ll . Most 
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comprehensive surveys of exact solutions to this nonlinear equation can be found in the handbook |13C 
it also describes a considerable number of generalized and functional separable solutions to nonUnear 
reaction-diffusion systems of two coupled equations without delay. 

The list of known exact solutions to equation ([Hi is quite limited. 

In general, equation ([T]) admits traveling-wave solutions, u = u{ax + /3t). Such solutions are dealt 
with in many studies (e.g., see the studies MiE and references in them). 

A complete group analysis of the non-linear differential-difference equation ([Jl was carried out 



in illh . Four equations of the form ([T} were found to admit invariant solutions; two of these equations 
are of limited interest, since they have degenerate solutions (linear in x). There was only one equation 
that involved an arbitrary function and had a non-degenerate solution: 

Ut = kuxx -'r u[—alnu + f{wu~ )], 6 = 6°"^, (2) 

where f{z) is an arbitrary function. The exact solution to this equation found in lilLI was 

u = exp(Cxe-''*)v?(t), (3) 

where C is an arbitrary constant and ip{t) is a function satisfying the delay ordinary differential equation 

if' it) = ^(t) [C'^ke-'"'' - aln ip{t) + /(^(t - T)ip'\t))] . (4) 



The other equation obtained in lllin that had a non-degenerate solution coincides, up to notation, with a 
special case of equations ^, at f{z) = ci + C2 In z. 

Remark 1 . It is noteworthy that equation (|2]i is explicitly dependent on the delay time r, which corresponds 
to a more general kinetic function F{u,w,t) than in equation ([TJ. It is only at a — in (l2]i that we have a 
kinetic function explicitly independent of r: F{u, w) — uf{w/u); in this case, © represents a separable solution, 

u = e^'^tp{t). 

In what follows, the term 'exact solution' with regard to nonlinear partial differential-difference 
equations, including delay partial differential equations, is used in the following cases: 

(i) the solution is expressible in terms of elementary functions or in closed form with definite or 
indefinite integrals; 

(ii) the solution is expressible in terms of solutions to ordinary differential or ordinary differential- 
difference equations (or systems of such equations); 

(iii) the solution is expressible in terms of solutions to linear partial differential equations. 

Combinations of cases (i)-(iii) are also allowed. 

This definition generalizes the notion of an exact solution used in JBOll with regard to nonlinear partial 
differential equations. 

Remark 2. Solution methods and various applications of linear and nonlinear ordinary differential-difference 
equations, which are much simpler than nonlinear partial differential-difference equations, can be found, for ex- 
ample, in 1.31.- 36.1. 

2 General description of the functional constraints method 

Consider a wide class of nonlinear delay reaction-diffusion equations: 

ut = kuxx + uf{z) + wg{z) + h{z), 

'W = u{x,t — T), Z = z{u,w), 

where f{z), g{z), and h{z) are arbitrary functions and z = z{u,w) is a given function. In addition, 
we will sometimes consider more complex equations where /, g, and h can additionally depend on the 
independent variables x or/and t explicitly. 



We look for generalized separable solutions of the form 

N 

U = J2^n{x)^n{t), (6) 

n=l 

where the functions ^n{x) and ^n{t) are to be determined in the analysis. 

Remark 3. For nonlinear partial differential equations, various modifications of the method of generalized 
separation of variables based on searching for solutions of the form (|6]l are detailed, for example, in ll27i BSlboll . 
These studies also present a large number of equations that admits generalized separable solutions. 

For nonlinear delay partial differential equations of the form (|5]l that involve arbitrary functions, the direct 
application of the method of generahzed separation of variables turns out to be ineffective. 

The new approach pursued in the present paper is based on searching for generalized separable 
solutions of the form ^ that satisfy one the following two additional /M«c?/o?ia/ constraints: 

z{u,w) = p{x), w = u(x,t — t); (7) 

z{u,w) = q{t), w = u{x,t — T). (8) 

These constraints represent difference equations in t with x playing the role of a free parameter. The 
function z = z{u, w) appears in equation dS]) as the argument of the arbitrary functions. The functions 
p{x) and q{t) are implicitly dependent on x and t (expressible in terms of ^n{x) and ^„(i), respec- 
tively) and are determined from the analysis of equations © or dSjl taking into account ©. It should be 
emphasized that there is no need to obtain general solutions of equations (|7]) or ([8]l; particular solutions 
will suffice. 

A solution to the difference equation (|7]) or ([8]), in view of Q, determines an allowed form of the 
exact solution, whose final representation is subsequently obtained by substituting the resulting solution 
into the specific equation Q. 

In what follows, constraints (Q and ^ will be referred to as afunctional constraint of the first kind 
and functional constraint of the second kind, respectively. 

For functional constraint of the first kind ©, the functions ^n(i) appearing in ^ are usually chosen 
in the form 

^n{t) = e^"* cos(/3„t), ^„(t) = e^"* sin(/3„t), 

with the parameters \n and pn determined from ^. 

Remark 4. We have introduced the term functional constraint by analogy with the term differential constraint, 
which is employed in the method of differential constraints used to seek exact solutions to nonlinear partial dif- 



ferential equations; for a description of this method and application examples, see, for example, ||3(XI37I - I39I1 ). In 
more complex cases, the original equation and functional constraints (jTjl and (|8]l can contain a time-varying delay 
T = T{t) instead of the constant delay t; see the comments after equation (|83] | in Section 8. 



Below we give examples of applying the above method to constructing generalized separable solu- 
tions to some equations of the form dS]) as well as more complex nonlinear partial differential equations 
with delay. 

3 The equation contains one arbitrary function dependent onw/u 

Equation 1. Consider the equation 

ut = kuxx + uf{w/u), (10) 

which is a special case of equation ^ with g = h = and z = w/u. 
1.1. The functional constraint of the second kind ([8]l becomes 

w/u = q{t), w = u{x,t — t). (11) 



It is clear that the difference equation (ITTT i can be satisfied with a simple separable solution 

u = ip{x)ip{t), (12) 

which gives q{t) = 'ijj{t — T)/ip{t). Substituting (fT2l ) into (ITOl ) and separating the variable, we arrive at 
the following equations for (p{x) and ^p{t)^. 

if" = aif, (13) 

v^'(t) = aki^it) + mim - r)/m), (14) 

where a is an arbitrary constant. 

The general solution of equation (IT3T l is expressed as 



(^(x) 



Ci cos(-y/|a[x) + C2 sin(-y/|a|x) if a < 0; 

Ci exp(— ^/ax) + C2exp(>/ox) if a > 0; (15) 

Cix + C2 if a = 0, 



where Ci and C2 are arbitrary constants. 

Remark 5. Simple separable solutions of the form (fT2] l are also admitted by the more general equation 

Uf = ku^x + uf{t,w/u), 

in which the kinetic function is additionally dependent on t. 

1.2. In the simple case p{x) = po = const, the functional constraint of the first kind (jT) for equation 
(ITOl) is written as 

w/u = po, w = u{x,t — t). (16) 

Solutions of the difference equation (IT6l ) with po > are sought in the form 

u = e^^vix, t), v{x, t) = v{x, t-r), (17) 

where c is an arbitrary constant and v{x, t) is a r-periodic function. In this case, w/u = po = e~'^^. 
Substituting dTT] ) into equation (ITOl ) yields a linear problem for determining v. 

vt = kvxx + hv, v{x,t) = v{x,t-T), (18) 

where b = /(e~^^) — c. 

The general solution to problem (ITST l. which will be denoted v = Vi{x,t;b) for convenience, is 
expressed as 

00 
Vi {x, t; 6) = ^ exp(-A„x) [An cos(/3„t - 7„2;) + S„ sin(/3„t - 7„x)] 

n=0 

00 

+ ^ exp(Artx) [Cn cos{/3nt + 7„x) + Dn sin(/3„t + 7„x)] , (19) 

n=l 

^" = -r' ^" = 1 — Tk — ) ' ^- = [ — 2k — ) ' ^^^^ 

where ^„, i3„, C„, and D„ are arbitrary constants such that the series ([19} with (120) as well as the 
derivatives {Vi)t and {Vi)xx are convergent; for example, the convergence can be ensured if one sets 
An = Bn = Cn = Dn = forn > A^, where N is some arbitrary positive integer. 

The following special cases can be distinguished: 

(i) r-periodic (in t) solutions to problem (fTSl ) decaying as x — >• oo are given by formulas (|T9] l 
and dlOll with Aq = Bq = 0, C„ = Dn = 0, and n = 1, 2, . . . ; 



(ii) r-periodic solutions to problem (ITST i bounded as x — )• oo are given by formulas (|T9] l and (|20] | 
with C„ = !)„ = and n = 1, 2, . . . ; 

(iii) a stationary solution is given by formulas (IT9l ) and (l20l) with ^d^, = Bn = Cn = Dn = and 
n=l, 2, .... 

To sum up, we have an exact solution to equation (ITOl ) 

7/ = e^Vi(x,t;6), 6 = /(e-^-)-c, (21) 

where c is an arbitrary constant and Vi{x, t;b) is a. r-periodic function determined by formulas ([19} and 
(Unil. 

1.3. Solutions of the difference equation (IT6] l with po < are sought in the form 

u = e'^^v{x,t), v{x,t) = —v{x,t — t), ill) 

where c is an arbitrary constant and v{x, t) is a r-antiperiodic function. In this case, w/u = po = —e~^'^. 
Substituting (|22] ) into (ITOl ) yields a linear problem for determining v: 

vt = kvxx + bv, v{x,t) = -v{x,t-T), (23) 

where b = f{—e~'^'^) — c. 

The general solution to problem (|23] ). which will be denoted v = V2{x,t;b) for convenience, is 
expressed as 

oo 

V2{x, t; 6) = ^ exp(-A„x) [An cos(/3„t - 7„j;) + Bn sin(/3„t - 7„j;)] 

n=l 

oo 

+ ^ exp(A„x) [C„ cos(/3„i + 7„x) + £)„ sin(/3„t + 7„x)] , (24) 

n=l 



2k J \ 2k 

where An, Bn, Cn, and Dn are arbitrary constants such that the series (l24l with (l25t as well as the 
derivatives {¥2)1 and {V2)xx are convergent. Solutions to problem (l23T l decaying as x — ;■ 00 and r- 
antiperiodic in t are given by formulas (l24l ) and (|25] ) with C„ = Dn = and n = 1, 2, . . . . 
To sum up, we have an exact solution to equation (ITOl) : 

u = e''^V2{x,t;b), 6 = /(-e"^^) - c, (26) 

where c is an arbitrary constant and V2{x, t; b) is a r-antiperiodic function determined by formulas (|24] | 
and(l25]). 



Remark 6. Solutions (fT9b-(l20l) and (I24l) -(l25b look very similar. However, the summation in the former solu- 
tion starts with n = rather than n = 1 and the expressions of /?„ are different. 

4 Equations with one arbitrary function dependent on a linear combina- 
tion of u and w 

Equation 2. Consider the equation 

ut = kuxx + bu + f{u-w), (27) 

which is a special case of equation ([5]) with f{z) = b,g = 0, and z = u — w (the function h has been 
renamed /). 



2.1. In this case, the functional constraint of the second kind ([8]l has the form 

u — w = q{t), w = u{x,t — T). (28) 

It is clear that the additive separable solution 

u = ip{x) + V(t) (29) 

satisfies the difference equation (1281 ). We have q{t) = 'ip{t) — 1^(1 — r). Substituting ( [291 ) into (|27] ) and 
separating the variables, one arrives at equations for determining (p{x) and 'ip{t): 

^Vxx + bip = a, (30) 

4{t) = bij{t) + a + /(V(t) - ^{t - r)), (31) 

where a is an arbitrary constant. 

The general solution to equation (l30l ) with 6 7^ and a = is expressed as 



I Ci cos(ax) + C2 sin(ax), a = \/b k if fe > 0; 

(^(x) = < (32) 

I Ci exp(— ax) + C2 exp(ax), a = yj—b/k if 6 < 0, 

where Ci and C2 are arbitrary constants. Solution (|29] l with 6 > is periodic in the space coordinate x. 
The general solution to equation (l30l) with 6 = and a / is 

v{x) = -^x^ + Cix + C2. (33) 

Remark 7. Solutions of the form(|29]l are also admitted by the more general equation 

ut = kuxx + hu + f{t, u - w), 
in which the kinetic function is additionally dependent on t. 

2.2. The functional constraint of the first kind © for equation (|27] ) has the form 

u — w = p{x), w = u{x,t — t). (34) 



The difference equation (1281 ) can be satisfied, for example, by choosing the generalized separable 
solution 

u = tip(x) + tpi^x), (35) 

which results in p{x) = Tip{x). 

Substituting (1351 ) into (|27]) yields ordinary differential equations for (^(x) and ip{x): 

kip'^, + b^ = 0, (36) 

fcV'L + bi^ + finp) - (^ = 0. (37) 

Equation (l36l ) coincides with (l30l ) at a = 0; its solution is given by formulas (|32] |. The nonhomogeneous 
constant-coefficient linear ordinary differential equation (|37] ) is easy to integrate. 

Using the above solutions (1291 ) and (1351 ) as well as the theorem below, one can obtain more complex 
exact solutions to equation (|27] ): these solutions can have any number of arbitrary parameters. 

Theorem 1 (nonlinear superposition of solutions). Suppose uq{x, t) is a solution to the nonlinear equa- 
tion ( 1271 ) and v = Vi (x, t; 6) is any t -periodic solution to the linear heat equation with source ( li8D . Then 
the sum 

u = UQ{x,t) + Vi{x,t]b) (38) 

is also a solution to equation ( 12 7D . The general form of the function Vi(x, t; b) is given by formulas ( Ii9l ) 
and^. 



Remark 8. For the nonlinear equation (|27] |. the particular solution uo{x, t) in formula (|38] | can be taken in the 
traveling-wave form uq — uo{ax + /3t). 

Remark 9. Formula (l38l l also remains valid for the more general equation 

ut = kuxx + hu + f{x, t,u- w), 

in which the kinetic function depends on three arguments. If / = f{x, u~w), then a stationary solution uq = uo{x) 
can be used as the first term in (l38l l. while the second term can be taken in the form (fT9Tl-(l20li. If / = /(t, u — w), 
then a spatially homogeneous solution uq — uo{t) can be taken for the first term in (l3Ft . 

Equation 3. Consider the equation 

Ut = kuxx + bu + f{u — aw), a > 0, (39) 

which is a special case of equation ^ with f{z) = b, g = 0, and z = u — aw; for convenience, the 
function h has been renamed /. 

3.1. The functional constraint of the first kind (|7]l for equation (1271 ) has the form 

u — aw=p{x), w = u{x,t — T). (40) 

The difference equation (l40l ) can be satisfied, for example, with the generalized separable solution 

u = e^^ipix) +ij{x), c = — Ina, (41) 

r 

which gives p{x) = (1 — a)ip{x). 

Substituting (1411 into (|39l ) leads to ordinary differential equations for ip{x) and tp{x): 

k^" + {b- c)(/7 = 0, (42) 

fcV" + ftV' + fiv) =0, ??=(!- a)V'- (43) 



Up to obvious renaming of variables, equation (I42l i coincides with (1301 ) at a = 0; its solution is given 
by formulas (|32] | where 6 must be substituted by 6 — c. 

3.2. Using the above solution (|4TI)-(|43]) as well as the theorem below, one can obtain more complex 
exact solutions to equation (l39l ): these solutions can have any number of arbitrary free parameters. 

Theorem 2 (generalization of Theorem [!}• Suppose uo{x, t) is a solution to the nonlinear equation f l?9l ) 
and V = Vi{x,t;b) is any T-periodic solution to the linear heat equation with source di^l ). Then the sum 

u = UQ{x,t) + e''^Vi{x,t;b- c), c=—\na, (44) 

r 



is also a solution to equation ([39\) . The general form of the function Vi{x, t; b) is given by formulas f lJPl 
and ( l20l ). 



Formula (I44l i makes it possible to obtain a wide class of exact solutions to the nonlinear equation ( 
Apart from (SB> one can take uq = uo{x) and uq = uo{t) as well as the more general, traveling wave 
solution no = 6{ax + /St) as the particular solution uq = uo{x, t); the constants a and /3 are arbitrary 
and the function 9{y) satisfies the ordinary differential-difference equation 

ka^e"{y) - (39'{y) + bd{y) + f{9{y) - ad{y - a)) = 0, y = ax + (3t, a = /3r. 
Remark 1 0. Formula ( l44l i also remains valid for the more general equation 

Ut = kuxx + bti + f{x, t,u ~ aw), 
in which the kinetic function depends on three arguments. 



Equation 4. Consider the equation 

ut = kuxx + bu + f{u + aw), a > 0, (45) 

which is a special case of equation ([5]) with f{z) = b,g = 0, and z = u + aw, for convenience, the 
function h has been renamed /. 

Theorem 3. Suppose uo{x,t) is a solution to the nonlinear equation d45D and v = V2{x,t;b) is any 
T-antiperiodic solution to the linear heat equation with source f l2JD . Then the sum 

u = uo{x,t) + e^V2{x,t;b - c), c=—lna, (46) 

r 



is also a solution to equation ([45\) . The general form of the function V2{x, t; b) is given by formulas (|24 
and fl25l). 



Formula (|46] | makes it possible to obtain a wide class of exact solutions to the nonlinear equation (l45T l. 
One can take uq = uo{x) and uq = uo{t) as well as the more general, traveling wave solution uq = 
9{ax + /3t) as the particular solution uo{x, t); the constants a and j3 are arbitrary and the function 9{y) 
satisfies the ordinary differential-difference equation 



r. 



ka^6"{y) - (39' {y) + b9{y) + f{9{y) + a9{y - a)) = 0, y = ax + (3t, a = (3 
Remark 1 1 . Formula ( |46] | also remains valid for the more general equation 

Ut = kuxx + bti + f{x, t,u + aw), 
in which the kinetic function depends on three arguments. 

5 Equations with two arbitrary functions dependent on a linear combi- 
nation of u and w 

Equation 5. Now consider the more general equation 

Ut = kuxx + uf{u — w) + wg{u ~ w) + h{u — w), (47) 

where f{z), g{z), and h{z) are arbitrary functions; in this case, either function f or g can be set equal to 
zero without loss of generality. 

5.1. The difference constraint © for equation Wf\ has the form (l34l ). The linear difference equa- 
tion (l34l ) can be satisfied, as previously, with a generalized separable solution of the form (l35T l. As a 
result, one can obtain equations for determining ip{x) and 4'{x)', these equations are not written out, 
since a significantly more general result will be presented below. 

5.2. The hnear difference equation (|34] | can be satisfied by setting 

'^ = y^[Pn{x) cos(/3„i) + Vn(x) sin(/3^t)] + t9{x) + i{x), Pn = ^^, (48) 

n=l 

where A^ is an arbitrary positive integer. In this case, we have p{x) = Tip{x) on the right-hand side of 
equation (l34l ). 

Substituting (1481 ) into WT\ and performing simple rearrangements, we obtain 

N 

Y,[^n cos(/3„t) + Bn sin(/3„t)] +Ct + D = Q, (49) 

n=l 

where the functional coefficients An, Bn, C, and D are dependent on ipn{x), ipn{x), 9{x), and ^(x) as 
well as their derivatives and independent of t. In (l49l ). equating all the functional coefficients with zero, 

8 



An = Bn = C = D = 0, we arrive at the following ordinary differential equation for the unknown 
functions: 

kif'; + ifnifire) + g{Te)] - pni^n = 0, 

fc< + MfirO) + gird)] + Pn^n = 0, 

ke" + e[f{Te) + g{Te)] = o, 

K" + if{rO) + {i- Te)g{Te) + h{T9) -6 = 0. 

Note that the third nonlinear equation admits the trivial particular solution 6 = 0; in this case, all 
other equations become linear with constant coefficients. 

Remark 1 2. Solutions of the form ( l48T l are also admitted by the more general equation 

Ut = kuxx + uf{x, M — w) + wg{x, u — w) + h{x, u — w). 
Equation 6. Consider the equation 

Ut = kuxx + uf{u — aw) + wg{u — aw) + h{u — aw), (50) 

where f{z), g{z), and h{z) are arbitrary functions, which generalizes equation (|39). 

6.1. The difference constraint © for equation (|50] | has the form (l40l ). The Unear difference equa- 
tion (l40l ) can be satisfied, as previously, with a generalized separable solution of the form (1411 . As a 
result, one arrives at equations for determining ip{x) and ip{x); these equations are not written here, 
since a much more general result is given below. 

6.2. The linear difference equation (|34] | can be satisfied by setting 



„ct 



1 



N 



+ X^['^n(a;)cos(/3„t) + ^„(x)sin(/3„t)] \ + ^{x 



n=l 



(51) 



c = — Ina, j3n 



27rn 



where N is an arbitrary positive integer In this case, we have p{x) = (1 — a)^(x) of the right-hand side 
of equation (l40l ). 

Substituting (ISTl ) into (l50l ) and reasoning in a similar fashion as for equation (|47] ). we arrive at the 
following equations for determining 9{x), ipn{x), ipn{x), and ^(x): 



A;0" + e 



1 



fiv) + -5(??) - c = 0, ?? = (!- a)i, 
L a J 

^<^n + 'fn fiv) + -g{v) - C 

L a J 

f{i) + -aiv) - c 

I a J 



ki'n + V'n 



- /^nV'n = 0, 
+ (3n^n = 0, 



Remark 1 3. Solutions of the form (fSTl i are also admitted by the more general equation 

Uf = kuxx + uf{x, u — aw) + wgi^x, u — aw) + h(x, u — aw), a > 0. 

Equation 7. Consider the equation 

Ut = kuxx + uf{u + aw) + wg{u + aw) + /i(n + aw), a > 0, 

where f{z), g{z), and h{z) are arbitrary functions, which generalizes equation (l45l . 
The difference constraint d?) for equation (|52] | has the form 



(52) 



u + aw = p{x). 



w = u{x,t — t). 



(53) 



The linear difference equation (l53T l can be satisfied by setting 



N 



^"^ ^["Priix) cos(/3„t) + iJnix) sin(/3„t)] + i{x), 



(54) 



n=l 

1 vr(2n + l) 

c = — m a, Pn = 

r r 



where A^ is an arbitrary positive integer. In this case, we have p{x) = (1 + a)^(x) on the right-hand side 
of equation (|53] ). 

Substituting (l54l) into (l52l ) and reasoning in a similar fashion as for equation (1471 ). we arrive at the 
following equations for determining ipn{x), ipn{x), and ^(x): 



1 



L a J 

ki'n + V'n /(f?) — g{v) - c 

L a J 



+ Pn^n = 0, 

ke + Cifiv) + 9{V)] + Kv) = 0, 7? = (1 + a)e. 

The last equation is independent. 

Remark 1 4. Solutions of the form ( l54l i are also admitted by the more general equation 

Uf — kuxx + uf(x, u + aw) + wg{x^ u + aw) + h{x, u + aio), a > 0. 

6 Equation with two arbitrary functions dependent on w^ + w' 

Equation 8. Now consider the equation 

ut = kuxx + uf{u^ + w^^) + wg{v? + w"^), 

which contains a nonlinear (quadratic) argument, z = u'^ + w'^. 
The difference constraint d?) for equation (l55T l has the form 

w^ + tii^ = p(x), w = u{x,t — t). 

The nonlinear difference equation (l53T l can be satisfied by setting 

u = ipn{x) cos(A„t) + il^nix) sin(A„i), 
TT(2n+l) 



(55) 



(56) 



(57) 



A 



2r 



n = 0, ±1, ±2, 



It is not difficult to verify that 

W = i-ir^nix) sin(Ant) + (-1)"+Vn(x) COs(Ant) 

and 

u^ + w^ = ipl{x) + ^^(x) = p{x). 

Substituting (|57] ) into (|55] ) and splitting the resulting expression with respect to cos(A„t) and sin(A„t), 
we arrive at a nonlinear system of ordinary differential equations for determining ipn{x) and tpn{x): 

kVn + ^nfivl + i^l) + (-1)"+Vn5(9'n + V'n) " AnV'n = 0, 
/C< + i^nf ivl + V'n) + (-l)>n5(V^n + V'') + AnV^n = 0. 

Remark 1 5. Solutions of the form (l54l l are also admitted by the more general equation 

Ut — kuxx + ^/[x, u + ui ) + wg{x, u + w ). 
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7 More general nonlinear delay partial differential equations 

Now consider nonlinear partial differential-difference equations of the more general form 

L[u]=M[u]+uf{z) + wg{z) + h{z), 

W = u(x,t — T), Z = z{u,w), 



where L and M are arbitrary constant-coefficient linear differential operators with respect to t and x: 

LM = i:^^^' MM=f;a.|^. (59) 



As before, f{z), g{z), and h{z) are arbitrary functions and z = z{u, w) is a given function. 

By setting L[u] = uu and M[u] = a?Uxx in (l58] ). we get the nonlinear delay Klein-Gordon equation 
(a delay hyperbolic PDE) 

utt = a^Uxx + uf{z) + wg{z) + h{z). (60) 

Remark 16. Oscillation properties of some hyperbolic partial differential equations with delay were studied, 
for example, in ll40l - l42ll . 



Many of the results obtained previously for nonlinear delay reaction-diffusion equations dSjl also 
apply to nonlinear delay partial differential equations of the general form (l58Tl-(l59Tl. 

Subsequently, in deriving determining equations, we use the following simple properties of the oper- 
ators L and M: 

L[(^(x)V(t)] = ^{x)l.[^{t)], L[(p(x) + V(t)] = L[^(t)], 
M[c^(x)V'(i)] = ^{t)M[^{x)], M[ip{x) + V(t)] = M[v9(x)]. 

Below we omit intermediate steps and give only final results. 

Equation 9. The equation 

l.[u]=M[u]+uf{w/u), (61) 

admits the multiplicative separable solution 

u = ip{x)ip{t), (62) 

with the functions 9?(x) and tpit) satisfying the linear ordinary differential and differential-difference 
equations 

M[<^] = op, (63) 

L[^(t)] = cm + i;{t)f{t, i,{t - T)/m) , (64) 

where c is an arbitrary constant. 

It should be noted that: 

(i) if c = 0, equation (|63] ) admits polynomial particular solutions; 

(ii) for any c, equation (|63] ) admits exponential particular solutions (p{x) = Ae^^, where A is an 
arbitrary constant and /3 is a root of the polynomial equation YliLi (^iP^ = c; 

(iii) equation (l64l ) admits exponential particular solutions ^(t) = i?e^*, where B is an arbitrary 
constant and A is a root of the transcendental equation 



5^M^ = c + /(e-^"). 



Equation 10. The equation 

L[u] = M[u\ + bu + f{u - w) (65) 
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admits additive separable solutions of the form 

M = V9(x) + V(t), (66) 

with the functions ip{x) and ^|^{t) described by the linear partial differential equation and ordinary 
differential-difference equation 

M[(p]=c- bip, (67) 

L[V;(t)] = c + bm + /(V^(t) - V(t - r)) (68) 

and c being an arbitrary constant. 

It should be noted that: 

(i) ifb = 0, equation (|67] ) admits polynomial particular solutions; 

(ii) if c = 0, equation (l67l) admits exponential particular solutions (p{x) = Ae^^, where A is an 
arbitrary constant and /3 is a root of the polynomial equation Y17=i '^if^^ ~ ~^' 

Using the above solution (l66l ) and the theorem below, one can find more complex exact solutions to 
equation (|65] ): these solutions can have any number of arbitrary free parameters. 



Theorem 4 (nonlinear superposition of solutions). Suppose uq{x, t) is a solution to the nonlinear equa- 
tion d65D and v = v{x, t) is any r-periodic solution to the linear equation 

\^[v\=M[v\+bv, v{x,t) =v{x,t-T). (69) 

Then the sum 

u = uo{x,t) + v{x,t) (70) 

is also solution to equation Ii65i . 

Apart from (l66l ). the particular solution uo{x,t) can be sought in the traveling-wave form uq = 
9{ax + fit). Particular solutions to the linear problem (|69) can be sought in the form 

N 

v{x, t) = ^ exp(-A„x) [An cos(/3„t - 7„a;) + Bn sin(/3„t - 7„x)] 

n=0 

N 

+ ^ exp(A„x) [Cn cos{(3nt + Tna:;) + Dn sin(/3„t + -/nx)] , (71) 

n=l 

Pn = 2lTn/T, 

where N is any positive integer, An, Bn, Cn, and Dn are arbitrary constants, and A„ and 7„ are constants 
determined from substituting (f/TI i into equation (1691 ). 

Equation 11. For the equation 

L[u] =M[u] + 6u + /(n-at(;), a > 0, (72) 

the following theorem holds. 

Theorem 5 (generalizes Theorem |4ll. Suppose uq{x, t) is a solution to the nonlinear equation ([72\) and 
V = v{x, t) is any T-periodic solution to the linear equation 

Li[v]=M[v] + bv, v{x,t) = v{x,t-T), (73) 

where Li[t;] = e^^*L[e'^*f] is a linear differential operator with constant coefficients. Then the sum 

u = uo{x,t) + e'^^v{x,t), c=— Ina, (74) 

T 



is also a solution to equation i\72}. 

Formula (1741 1 makes it possible to obtain a wide class of exact solutions to the nonlinear equation d72] i. 
Simple solutions of the form uq = uq{x) and uq = no(t) as well as the more general, traveling wave 
solution Uq = 6{ax + fit) can be taken as the particular solution UQ{x,t). Just as before, particular 
solutions to the linear problem (TTBT i can be sought in the form (TtTI ). 
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Equation 12. For the equation 

h[u]=M{u] + hu + f{u + aw), a > 0, (75) 

the following theorem holds. 

Theorem 6. Suppose UQ{x,t) is a solution to the nonlinear equation d75D and v = v{x,t) is any r- 
antiperiodic solution to the linear equation 

'Li[v\=M[v\ + bv, v{x,t) = -v{x,t-T), (76) 

where Li [v] = e^'^^L[ef^v] is a linear differential operator with constant coefficients. Then the sum 

u = UQ{x,t) + e'^*'v{x,t), c=— Ina, (77) 

r 

is also a solution to equation d75D . 

Formula (TTTT i makes it possible to obtain a wide class of exact solutions to the nonlinear equation 
(IT?] ). Simple solutions of the form uq = uq{x) and uq = uo{t) as well as the more general, traveling 
wave solution uq = 9{ax + /3t) can be taken as the particular solution uo{x, t). Particular solutions to 
the Unear problem (1761 ) can be sought in the form 

N 

v{x, t) = ^ exp(-A„x) [An cos(/3„t - 7„a;) + 5„ sin(/3„t - 7„x)] 

n=l 

N 

+ ^ exp(A„x) [Cn cos{(3nt + Jnx) + Dn sin(/3„t + 7„x)] , (78) 

71=1 

7r(2n-l) 

Pn = 



r 



where A^ is any positive integer, An, Bn, Cn, and Dn are arbitrary constants, and A„ and 7„ are constants 
determined from substituting (iTST i into (TTSl l. 

Equation 13. Exact solutions to the nonlinear delay partial differential equation 

L[n] = M[u] + uf{x, u — w) + wg{x, u — w) + h{x, u — w) (79) 

are sought in the form (@8). 

Equation 14. Exact solutions to the nonlinear delay partial differential equation 

'L[u\=y[[u\ + uf{x,u — aw) + wg{x,u — aw) + h{x,u — aw), a > 0, (80) 

are sought in the form (l5Tb . 

Equation 15. Exact solutions to the nonlinear delay partial differential equation 

L[ti] = M[ti] + uf{x, u + at(;) + wg{x, u + aw;) + h{x, u + atf), a > 0, (81) 

are sought in the form (|54). 

Equation 16. Exact solutions to the nonlinear delay partial differential equation 

L[n] = M[u\ + u/(x, -u^ + w"^) + 'u;5'(x, -u^ + w'^) (82) 

are sought in the form (l57b . 
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Remark 1 7. In equations < [79l t, dSOl t, dSTT t, and <l82b , the linear operator M can depend on x and have the form 

ni Q^ 

i—l 

where ai{x) are arbitrary functions. 

Remark 18. In equations ( l79b . dSOb , dSTT l. and ( |82] |. the functions /, g, and h can expUcitly depend on several 
space coordinates, x = {xi, . . . ,Xs), and M can be an arbitrary linear differential operator of any order (first, 
second, or higher) with respect to the space coordinates whose coefficients can be explicitly dependent on a; = 
{xi, . . . , Xs)- Just as before, exact solutions to such equations are sought in the form (|48] |. (fSTT l. (l54l i. and (|57] |. 
respectively, where x is understood as a; = {xi, . . . , Xg)- 

In particular, one can choose an elliptic operator. 



which arises in the studies 1I20II21I1 in the case q = 0. 



"^ dxi ''^ dxj ^-^ * dxi 

zj — 1 ■' z— 1 



8 General partial differential equations with time-varying delay 

In very much the same manner, one can obtain exact solutions to some nonlinear partial functional 
differential equations with time- varying delay of general form. 

Below we give a few simple examples to illustrate the aforesaid. 

Equation 17. Consider nonlinear partial functional differential equations of the form 

l.[u]=M[u]+uf{w/u), w = u{xX{t)), (83) 



where C,{t) is a given function. In applications (e.g., see 1201] '). it is conventional to write C,{t) = t — T{t), 
where T(t) is a time-varying delay such that < T{t) < tq. 
In this case, the functional constraint dS]) becomes 

w/u = q{t), w = u{x, C(t)). (84) 

The functional equation (l84l ) can be satisfied with the simple multiplicative separable solution 

u = ^{x)i^{t), (85) 

which gives q{t) = 'ip{({t))/ip{t). Substituting (|85] ) into (|83] ) and separating the variables, one arrives at 
the following equations for (/?(x) and Tp{t): 

M[^] = ap, 

L[V'(t)] = cV^(t) + V(t)/(i,V'(C)Mt)), C = C{t), 

where c is an arbitrary constant. 

Equation 18. Consider nonlinear partial functional differential equations of the form 

L[u]=M[u] + bu + f{u-w), w = u{xX{t))- (86) 

In this case, the functional constraint dS]) becomes 

u — w = q{t), w = ui^xXit))- (87) 

The equation admits an additive separable solution of the form 

u = ip{x) + ip{t), (88) 

which gives q{t) = ifj^t) — ip{C,{t)). Substituting (|88] ) into (l86l ) and separating the variables, one arrives 
at the following equations for ip{x) and ipit): 

M[{p\ = c-b(f, 

L^)] = c + 6^(t) + fim - HO), c = m, 

where c is an arbitrary constant. 
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9 Brief conclusions 

To summarize, we have proposed a new method, which we call the functional constraints method, for 
constructing exact solutions of nonlinear delay reaction-diffusion equations of the form 

ut = kuxx + F{u,w), 

where u = u{x, t), w = u{x, t — r), and r is the relaxation time. The method is based on searching for 
generalized separable solutions of the form 



TV 
U = 

n=l 



J2^n(.x)rin{t), 



with the functions ^n(x) and r/„(t) determined from additional functional constraints (difference or func- 
tional equations) and the original delay partial differential equation. 

All of the equations considered contain one or two arbitrary functions of a single argument. We have 
described a considerable number of new exact generalized separable solutions and a few more complex 
solutions representing a nonlinear superposition of generalized separable solutions and traveling wave 
solutions. 

All solutions involve free parameters (in some cases, infinitely many parameters) and so can be suit- 
able for solving certain problems and testing approximate analytical and numerical methods for nonlinear 
delay PDEs. 

The results have been extended to the following classes of equations: 

(i) nonlinear delay partial differential equations of the form 

L[u] =M[u] +F{u,w), 

where L and M are arbitrary constant-coefficient linear differential operators of any order with respect 
to the independent variables x and t; in particular, this broad class includes the nonlinear delay Klein- 
Gordon equation, which is a hyperbolic-type delay equation; 

(ii) some nonlinear delay partial differential equations that are explicitly dependent on either of the 
independent variables, x or t; and 

(iii) some nonlinear partial functional differential equations with time- varying delay of general form. 
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